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Identifying the source parameters from a gravitational-wave measurement alone is limited by our ability to 
discriminate signals from different sources and the accuracy of the waveform family employed in the search. 
Here we address both issues in the framework of an adapted coordinate system that allows for linear Fisher- 
matrix type calculations of waveform differences that are both accurate and computationally very efficient. We 
investigate statistical errors by using principal component analysis of the post-Newtonian (PN) expansion co- 
efficients, which is well conditioned despite the Fisher matrix becoming ill-conditioned for larger numbers 
of parameters. We identify which combinations of physical parameters are most effectively measured by 
gravitational-wave detectors for systems of neutron stars and black holes with aligned spin. We confirm the 
expectation that the dominant parameter of the inspiral waveform is the chirp mass. The next dominant parame- 
ter depends on a combination of the spin and the symmetric mass ratio. In addition, we can study the systematic 
effect of various spin contributions to the PN phasing within the same parametrization, showing that the inclu- 
sion of spin-orbit corrections up to next-to-leading order, but not necessarily of spin-spin contributions, is crucial 
for an accurate inspiral waveform model. This understanding of the waveform structure throughout the param- 
eter space is important to set up an efficient search strategy and correctly interpret future gravitational-wave 
observations. 

PACS numbers: 04.8().Nn, ()4.25.Nx, 04.30.Db, 
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I. INTRODUCTION 

Ground-based gravitational-wave (GW) detectors of 
the Laser Interferometer Gravitational-wave Observatory 
(LIGO) [1-4] and Virgo [5, 6] collaborations are currently 
being upgraded to provide sensitivities capable of directly 
detecting GWs from compact binary coalescences of binary 
black holes (BHs) and neutron star (NS) systems [7]. Such 
detections would constitute the first direct detection of 
NS-BH and binary BH systems. The gravitational waveforms 
from these systems will provide unprecedented information 
about the physical nature of these systems, and extracting 
these information relies on overlapping the noisy detector 
data with accurate theoretical signal predictions. 

The waveform from an inspiraling compact binary system 
can be obtained from knowledge of the energy and energy 
flux of the system. In general relativity these can be calcu- 
lated perturbatively in a v/c expansion (where v is the relative 
velocity of the bodies, c is the speed of light), known as a 
post-Newtonian (PN) expansion (see, e.g., [8] and references 
therein). These calculations provide the coefficients in such 
an expansion in terms of the fundamental physical parameters. 
Various different expansion schemes exist that lead to differ- 
ent approximants [9]. For quasi-circular, adiabatic orbits, the 
tangential velocity v can be related to the orbital frequency 
of the compact bodies, which, for the dominant GW mode, is 
equivalent to half the GW frequency. In this way an expansion 
in GW frequency, /, can be obtained, with each successive 
term corresponding to a higher order in the PN expansion. 

These waveforms depend on a number of physical parame- 
ters such as the masses and magnitudes and orientations of the 
objects' spins. An important task will be extracting as much 
of this information as possible given the observational con- 



straints and detector sensitivities. Although the masses and 
spins of the constituent objects are typically the parameters of 
greatest astrophysical interest, in practice the detectors are ac- 
tually sensitive to combinations of these parameters. This is 
because only sufficiently strong waveform variations, rather 
than variations in the source parameters, can be distinguished 
by the detector An example of this may be seen already in the 
Newtonian regime where the waveform of the binary inspi- 
ral depends only on what is commonly called the chirp mass, 
M., a combination of the two individual masses mi and m2, 
given by Al =: {71111712)'^^^ /{mi +7712)^^^. Systems with the 
same chirp mass emit the same GW phase in this simple ap- 
proximation, and they cannot be distinguished. Although this 
degeneracy does not hold in general, it remains true that GW 
detections can put much stronger constraints on a combina- 
tion of the masses characteristic of the binary system, in this 
case the chirp mass, than it can on the physical parameters of 
either individual object. 

One major difference between Newtonian dynamics and 
general relativity is that in the relativistic domain the spin an- 
gular momenta of the inspiraling objects affects the orbit and 
thus the gravitational waveform. In general relativity these 
spin effects first show up in the PN expansion as spin-orbit in- 
teractions at 1.5PN order, corresponding to (w/c)^ order cor- 
rections. At higher orders one also encounters spin-spin inter- 
actions. With PN corrections beyond the Newtonian term and 
in particular spin effects, it is less obvious which combinations 
of the physical parameters are most accurately measured, but 
empirical overlap studies [10, 11] have recently pointed out 
that apart from the chirp mass, there is a close degeneracy be- 
tween mass ratio and spins for BHs with spins aligned with 
the orbital angular momentum. (The waveforms of spinning 
NSs exhibit additional degeneracies, e.g., between the NS spin 



and the equation-of-state dependent quadrupole moment, but 
in this paper we neglect the spin of the NS as it is expected to 
be small in compact binary systems [12, 13].) 

Here we formalize the search for well-measurable parame- 
ters and degeneracies in the PN inspiral waveform. We em- 
ploy a linearization of waveform differences equivalent to 
the Fisher-matrix approximation [14], but we demonstrate 
that a convenient higher-dimensional coordinate choice in 
terms of the PN expansion coefficients allows for very accu- 
rate, yet computationally cheap calculations of the waveform 
(dis)agreement. 

The method we employ is similar to the one used by 
Tagoshi and Tanaka [15], Sathyaprakash and Schutz [16], Pai 
and Arun [17] and Brown et al. [18]. We write the waveform 
as a series expansion in frequency space and use the expansion 
coefficients as model parameters to construct a Fisher matrix. 
Using the eigenvalues and eigenvectors of this Fisher matrix 
we then determine which combinations of the expansion co- 
efficients the detector is most sensitive to, which amounts to 
finding the principal components of the Fisher matrix. In con- 
trast to Pai and Arun our focus is determining the best mea- 
sured combinations of parameters given aligned spinning gen- 
eral relativity waveforms and an Advanced LIGO noise curve 
[19]. We also discuss implications of a parameter dependent 
frequency cutoff. 

An example of principal component analysis (PCA) of 
spinning signals for the proposed LISA (Laser Interferometer 
Space Antenna) space-based mission was given in [16]. We, 
however, consider ground-based detectors, specifically Ad- 
vanced LIGO, where the expected signal-to-noise ratio (SNR) 
of most detections is rather low. Then it will be especially 
important to know which parameters can be measured since 
it is unlikely that all physical parameters will be measurable 
with reasonable accuracy. The number of principal compo- 
nents with a prescribed accuracy determined by the PCA will 
define an effective dimension of the space of solutions to be 
searched [18]. A solution space with a small effective dimen- 
sion will need relatively few templates to be searched, which 
speeds up search times [15, 18]. 

Our main aim here is to determine what the principal com- 
ponents represent physically. This cannot reduce any uncer- 
tainty in the measurement of physical parameters, which is 
typically large because of the correlation between the param- 
eters. However, an uncorrected set of parameters will give 
more tightly constrained directions in the likelihood space and 
also provide a convenient coordinate system in which to eval- 
uate the overlap of differing waveforms. 

In calculating the principal components we use as much in- 
formation about the PN coefficients of aligned spinning sys- 
tems as is currently available. The functional dependence 
of the PN coefficients on the physical parameters dictates 
how our principal components vary across parameter space. 
Furthermore, we investigate what the contribution of vari- 
ous terms in these coefficients are and how excluding them 
might affect parameter values through parameter bias. This 
helps to show which terms are important in the parameter es- 
timation problem and gives some indication of how yet-to-be- 
calculated terms may affect our results. 



This paper is organized as follows. After a general in- 
troduction of the Fisher-matrix approximation and PCA in 
Sec. II A, we specify the waveform model in Sec. II B and ar- 
gue by virtue of the Bauer- Fike theorem that even though the 
higher-dimensional Fisher matrix may be ill-conditioned, the 
corresponding principal components with large eigenvalues 
can be calculated accurately. In Sec. Ill we demonstrate ex- 
plicitly the superiority of our approach over standard Fisher- 
matrix estimates in terms of physical parameters, and we ex- 
tensively analyze the physical dependence of the first principal 
components. Sec. IV extends our algorithm to the case of dif- 
ferent waveform models, which enables us to identify crucial 
contributions to the PN phasing. We conclude with Sec. V. 



II. WAVEFORM AND METHODOLOGY 

A. Fisher matrix and principal component analysis 

The fundamental question that is underlying matched-filter 
searches for GWs is how different is a waveform hi (the de- 
tected signal) to another waveform ft,2 (the template). Assum- 
ing a GW detector with noise spectral density S'„(/), the ap- 
propriate difference is commonly defined by 



11^1 - ^2JP = {hi - /l2, hi - /l2), 



where the inner product reads 



(,,,.,,. 4Ro£"Wfl,,^ 



(1) 



(2) 



Here, h denotes the Fourier transform of h and * is the com- 
plex conjugation. Throughout this paper, we shall assume the 
noise spectral density of Advanced LIGO, in the zero detuned 
high power configuration detailed in [19], with a lower cut- 
off frequency at /mjn ==15 Hz. The upper frequency cutoff, 
.fmax^ is determined by the signal templates we assume. Since 
we are dealing exclusively with inspiral signals and ignoring 
the merger and ringdown phase, it is common practice to cut 
the inspiral frequencies at the equivalent innermost stable cir- 
cular orbit of the Schwarzschild spacetime, at 



/n 



l/(63/2^M) 



(3) 



(where M is the total mass of the system). We could use a 
more general, possibly spin-dependent, description of the up- 
per cutoff frequency, but that is not the focus of this paper, and 
our algorithm is readily applicable to more complicated forms 

or ./max- 

By evaluating the distance (1) or the overlap (2) one can, 
with some confidence, draw conclusions about the origin of 
the detected signal, if the template agrees very well with the 
data. Two main issues arise, however. Various templates with 
different source parameters may all agree well in terms of the 
predicted GW signals, in which case it would be impossible 
to definitively identify the true source parameters from just 
the GW observation. In addition, the template family may not 
be an exact representation of the real waveform, which again 



limits our ability to unambiguously identify the source of the 
detected signal. 

We shall analyze both effects here, statistical uncertainties 
due to similar waveforms for different parameters and system- 
atic biases due to inaccuracies in the waveform model. As 
such analyses are crucial for the correct interpretation of GW 
observations, there are already a number of publications ad- 
dressing these issues under various assumptions. In particular, 
PN inspiral waveforms have been analyzed by several authors, 
who calculated the measurement accuracy of various source 
parameters assuming a similar frequency -domain model to the 
one we employ here, but to lower PN expansion order [20-24] 
or only for nonspinning systems [25]. Systematic errors be- 
tween different PN approximants describing nonspinning sys- 
tems have been studied with extensive overlap calculations in 
[9]. 

We elaborate on the existing insights here by improving 
the linearization of (1) through a suitable coordinate choice 
in combination with a PCA, which allows us to understand in 
a systematic way which combinations of physical parameters 
are best constrained and which analytical contributions to the 
inspiral waveform are crucial to correctly recover the source 
parameters. 

The details of our strategy are as follows. We assume that 
hi and ft,2 can be parametrized by a common waveform man- 
ifold h such that hi = h{9) and ft,2 = h{6), where 9 is the 
vector of waveform parameters with components 6i (we shall 
put these in concrete terms in Sec. II B). With a minimization 
of the distance (1) in mind, we next apply the well-known lin- 
earization 



\h{e)-h{6)r^Y.^. 



jj A6'i A0j , 



(4) 



«j 



where A0 = 9 — 6 and F,, is the Fisher information matrix. 



_ / dh^ dh 



(5) 



describes the principal components of the system. 

Working with these coordinates rather than the original 
parametrization has the great advantage that the waveform dif- 
ference (4) becomes simply 



\hie)-hi9)r = J2K^i^l 



(8) 



We can now easily conclude from the size of the eigenvalues 
which principal components (or which principal directions in 
parameter space) affect the waveform strongly. This will be 
important to understand how well constrained and therefore 
measurable certain parameter combinations are, given that 
waveforms that differ below the noise floor cannot be distin- 
guished from each other. 

Typically, the smallest difference that is measurable is 
quoted to be 



\h{9)-hi9)\\ 



1, 



(9) 



see [26] and further discussions in [27, 28]. Here, however, we 
follow the recent discussion by Baird et al. [10], who detail 
that the 90% confidence interval in the posterior probability 
distribution is given to linear order by 



\h{9) - hi9)f < xl 



(10) 



where xl is the x^ value for which the probability of obtaining 
that value or less in a x^ -distribution with k degrees of free- 
dom is 90%. We shall later restrict ourselves to three physical 
parameters (the two masses of the compact objects and one 
spin magnitude) where x§ « 6.25 and waveforms with dis- 
tance 



\h{9)-h{9)f <6.25 



(11) 



cannot be distinguished at 90% confidence. Note that for 
SNR 10, this is approximately equivalent to the region of 
waveforms with overlap greater than 97% [10]. 



For more details about this approach and the validity of the 
Fisher-matrix approximation, see for instance [14]. We sim- 
ply use it as a convenient linearization here. 

The inverse of the Fisher matrix is the covariance matrix of 
the waveform parameters, and instead of quoting the variances 
of the used parameters (as done, e.g., in [20-25]), we proceed 
by diagonalizing the Fisher matrix. The result can be written 
as 



r,: 



22 ^fk ^kl Mj, 
k,l 



(6) 



where Aij denotes the j-th component of the z-th eigenvector 
of the Fisher matrix. S^j is a diagonal matrix with positive 
eigenvalues on the diagonal, i.e., E^j = XiSij. Since the 
eigenvectors are also eigenvectors of the covariance matrix, 
we have thus performed a PCA, and the vector 



Mi 



Ea. 



(7) 



B. Waveform model 

Let us specify in this section what waveform manifold h{9) 
we are using, which set of parameters 9i we are considering, 
and how we can take advantage of the methods presented in 
the previous section. 

In GW data analysis, the inspiral waveform of a coalesc- 
ing compact binary is most conveniently expressed in a closed 
form in the Fourier domain, which allows for a very fast eval- 
uation of thousands of templates, if needed. We shall use the 
same signal description here, which is commonly referred to 
as the "TaylorF2" approximant. Derived from the stationary- 
phase-approximation of the PN energy balance law [29, 30], 
it reads 



h{J)=Ar^''^e'^'^f\ 



(12) 



where A is an amplitude term which we set by requiring a par- 
ticular SNR, (ft,, h) = p^. We do not consider contributions 



from higher harmonics, which can be found in [31, 32]; hence 
A is simply a constant determined by the binary's total mass, 
distance, orientation and sky location. With this assumption 
(12) is often called the restricted PN waveform [9]. 

The phase, ip(f), is given as a series in the gravitational 
wave frequency /, 



Hf) 



N 

E 

fc=0 



(fe-5)/3 ^ 



■^fe + ^i°s log(///o) 



(13) 



where /g is an arbitrary reference frequency we introduce 
to make all coefficients dimensionless. The expansion co- 
efficients t/jk and ip^^ have been determined within the PN 
formalism in various publications (see [8] and references 
therein). Currently, the highest PN order to which they are 
known is 3.5 (k — 7) for nonspinning contributions. Spin- 
dependent contributions enter as leading-order and next-to- 
leading order spin-orbit effects at 1.5 and 2.5PN order (k = 
3, 5) [33-35]. We also include the tail-induced spin-orbit con- 
tribution to the flux at 3PN order (fc == 6) [36]. In addition, 
spin-spin effects are included at relative 2PN order [33, 37- 
39] and a 2.5PN contribution to the flux is associated with 
energy flux into the BH [40]. A collection of explicit expres- 
sions for the PN phase coefficients we employ can be found 
either in Eq. (A21) of [41] or Eq. (2.91) of [42], in combina- 
tion with Eq. (6.6) of [36]. 

Three physical source parameters define the binaries we 
consider: one object is a BH of mass mi, the other is a NS 
of mass ?Ti2. In addition, the BH is allowed to have a spin Si 
aligned with the direction of the orbital angular momentum L 
of the binary, which we parametrize through the dimension- 
less quantity 



Xi 



L-Si 

^2 



(14) 



The NS spin is assumed to be negligible and set to zero. We 
could also include the spin of the second compact object with- 
out any modification to our algorithm, as long as the spins are 
aligned with the orbital angular momentum. However, astro- 
physical expectations are that NSs in compact binaries do not 
spin rapidly [12, 13], which is confirmed by the fact that the 
highest NS spin parameter in a binary observed to date [43] 
has a value of x ^ 0.05. In addition, it was argued from a 
purely GW data analysis point of view that two spin parame- 
ters can efficiently be mapped onto one effective spin parame- 
ter, essentially without changing the waveform manifold [44- 
46]. In that sense, one spinning BH can simply be seen as 
a representative of the class of aligned-spin systems with the 
same effective spin. 

We thus consider the PN coefficients in (13) generally as 
functions of {mi, m2, xi)- Note, however, that we work un- 
der the assumption of general relativity, in which ipi and many 
log terms are exactly zero.' 



Two coefficients deserve further attention, ip^ is a constant 
phase that comes with no frequency-dependent factor, and it 
represents what is usually referred to as an additional param- 
eter: the "phase at coalescence" or the initial phase. The cor- 
responding initial time or "time at coalescence" is included in 
V's, which is a linear contribution to the phase. When we esti- 
mate waveform differences later, we are not interested in any 
discrepancy caused by time or phase shifts, so we shall project 
these parameters out of the Fisher matrix. 

Now that we have defined our parameters and waveform 
model, we could proceed by calculating the Fisher matrix 
(5) for the parameters {mi, m,2,Xii 4>o, to), where (f>o and to 
represent the free time and phase shift, or equivalently for 
any other 5-dimensional parametrization that can uniquely be 
mapped to the above parameters. It was found, however, in 
various publications that assuming a moderate SNR (between 
10 and 20) and the noise spectral density of the (initial or 
advanced) LIGO detector leads to rather large statistical pa- 
rameter biases in some parameters of interest. The lineariza- 
tion of the waveform difference implies that waveforms can 
definitively be told apart only if parameters like the individual 
masses or the mass ratio are considerably different [22-24] 
(to the order of tens of percents and more). The validity of 
the linearization is rightly questioned in these cases, and the 
true confidence interval (10) likely has more structure than the 
ellipsoid predicted in the Fisher-matrix approach. 

We partly circumvent these problems here by using the PN 
coefficients themselves as free parameters rather than the orig- 
inal physical parameters, i.e.. 



0.^{^u}u{i.'°^}. 



(15) 



All nonzero PN coefficients up to iV — 8 are now consid- 
ered to be free parameters, which yields a ten-dimensional 
parameter space that includes eight tpk (k = 0, 2, 3, ... , 8) 
and two tp^'^ (k = 5, 6). As noted above, the same idea, 
combined with an eigenvector analysis, has already been pre- 
sented by Tagoshi and Tanaka [15] and Brown et al. [18] in 
the context of template placing, and by Pai and Arun [17] to 
probe PN theory with GW observations. We shall show be- 
low how this trick is also useful to understand statistical and 
systematic errors of inspiral waveforms. One important fea- 
ture we will exploit is that the Fisher matrix (5) becomes al- 
most parameter-independent for the choice (15), which makes 
extrapolating the waveform difference (4) to large parameter 
variations much more accurate [15, 49]. Specifically, the only 
parameter dependence in the resulting matrix 



\A\' 



log'(///o) 

Sn{f) 



df 



(16) 



' In non-Einstein theories some of tliese terms may not be zero, for example 



in Brans-Dicke gravity a term arises already at a value of fc = —2 [47]. 
In order to capture waveforms of any gravitational theory, as performed 
in for example [48], who include a possible non-zero k = 1 term, these 
terms would need to be incorporated. In this work we focus on the general 
relativity results and set these terms explicitly to zero. 
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FIG. 1. Geometric interpretation of calculating waveform differences 
in a higher-dimensional flat space (8). The actual manifold of inspi- 
ral waveforms is illustrated as a two-dimensional curved surface that 
can be embedded in a higher-dimensional flat space [see coordinate 
choice (15)], here depicted in three dimensions. The confidence in- 
terval has a trivial geometry in the higher-dimensional space (here 
illustrated as a ball around the target waveform, which is shown as a 
white dot), which has to be projected back onto the physical wave- 
form manifold. 



value of Tij, can become quite large for the matrices we are 
considering. Calculating the eigenvalues for these matrices, 
however, is still a well-conditioned problem by virtue of the 
Bauer-Fike theorem which states [51] 



|Ar-Ar+5r|<«2(A)pr|l, 



(18) 



for a diagonalizable matrix F with eigenvalue Ar and a per- 
turbed matrix F + (5F with eigenvalue Ar+ar- The factor K2 
is the condition number of the diagonalizing matrix A and 
||(5F||2 denotes the matrix 2-norm of the perturbation matrix 
(5F. Since in the case of a real symmetric matrix F, the di- 
agonalizing matrix A is orthogonal with condition number 
K2(A) = 1 we are guaranteed that the eigenvalue problem 
is well-conditioned. A small change in F caused by numerical 
error will only cause a small change in the eigenvalues. In fact 
this is true for any normal matrix A satisfying A^A — AA^ 
with well-separated eigenvalues. This means that the large 
eigenvalues (and this is what we are interested in) of the Fisher 
matrix can be reliably computed even when the original ma- 
trix contains many PN coefficients and is itself ill-conditioned. 
Furthermore, under a perturbation of the Fisher matrix by 
(5F, the eigenvectors, Xi, are perturbed by Sxi, given in our 
case by [51] 



is inherited from the upper cutoff frequency (3). The expo- 
nents K and S are solely functions of i and j. 

Through this convenient parametrization, finding the con- 
fidence interval around a given target signal transforms to a 
simple geometric exercise in flat space. An illustration of ad- 
vantages and caveats of this change to higher-dimensional co- 
ordinates is given in Fig. 1. 

As mentioned earlier, we need to project the time and phase 
shifts out of the Fisher matrix, and we can do so following the 
simple procedure [50] 



r,: 



a.b 



r„7"'r, 



bj- 



(17) 



Here, 7"^ denotes the inverse of the two-dimensional subma- 
trix of F coiTesponding to the parameters ijj^ (phase shift) and 
-08 (time shift). The projected matrix Tij is effectively eight- 
dimensional, thus we shall find eight eigenvectors and eigen- 
values that govern waveform changes according to (8). 

Let us mention two caveats before we continue with pre- 
senting our results. First, in calculating the Fisher matrix with 
more than just our five physical parameters, we implicitly treat 
all PN coefficients as independent, which is clearly not true for 
the underlying lower-dimensional waveform manifold. So we 
have to take care to "project" our results back onto physically 
meaningful quantities. We shall do so by re-expressing the 
principal components not just as functions of PN coefficients, 
but eventually as functions of the physical parameters. 



6xi = ^ 



j¥'i 



xj STxi 
(A. -A,)- 



o{\\sn 



(19) 



For the eigenvectors with large eigenvalues A^, this ensures 
that their components are also well-conditioned despite the 
matrix F being ill-conditioned. This is sufficient for our pur- 
poses since we focus on waveform differences calculated via 
(8), which are dominated by the eigenvectors with large eigen- 
values. In numerical testing it was observed that the large 
eigenvalues are always well separated and that the eigenvalues 
and eigenvector components were stable for small changes of 
the Fisher matrix F. This stability of the diagonalization pro- 
cess applies to both the numerical eiTors and also to small 
deviations in the noise spectral density, which will also give 
rise to small changes in the Fisher matrix. 



III. MEASURABLE PARAMETER COMBINATIONS 
STATISTICAL ERROR 



Let us now present the results we obtain with the method 
and waveform model introduced in Sec. II. We first estimate 
statistical errors by asking the question which waveforms in a 
neighborhood of a given signal are similar to the latter to an 
extent that they cannot be distinguished at a 90% confidence 
level? We shall show that 



C. Numerical stability 

The other issue to be aware of is that the condition number, 
which we can view as the ratio of greatest and smallest eigen- 



(1) our method of estimating waveform differences is supe- 
rior to standard Fisher-matrix estimates (that are carried 
out in terms of the physical parameters), and we find no 
considerable differences between our computationally 
cheap method and full overlap calculations; 



(2) because of an approximate degeneracy between mass 
ratio and spin, the individual masses cannot be deter- 
mined accurately in GW observations alone, but 

(3) through a PCA we are able to identify the parameter 
combinations that are accurately measurable. 

Our results complement previous publications that estimated 
the measurability of source parameters of spinning systems 
either by Fisher-matrix calculations [22-24] or recently by di- 
rect overlap calculations [10, 11]. 



A. Advantage of different parametrizations 

The target system we consider for illustration is a binary 
containing a nonspinning NS of mass 1.35Mq and a BH with 
5Mq and a spin of xi — 0.3. We further assume a moder- 
ately high SNR of 20. We can now easily demonstrate the 
efficacy of our approach by estimating the 90% confidence in- 
terval [defined by (10)] around the target signal with various 
strategies and compare the results with proper overlap calcu- 
lations. 

The standard Fisher-matrix estimate in terms of 9i = 
{\ogM,logri,xi,ta,(j)o} is the simplest way of obtaining a 
multidimensional ellipse around the target parameters. Here 
we adopt the commonly used parametrization in terms of the 
symmetric mass ratio rj and the total mass M, 



mass ratio 
5 4 3 



V = 



mi 1712 



M = mi + 7712 



(20) 



Using the logarithms of the mass-dependent parameters im- 
proves the condition number of the Fisher matrix, and the 
square root of the diagonal elements of the inverse matrix 
(properly scaled) immediately yield the dimensions of the 
confidence interval. 



100%, 



AM ^ Ari 

— — « 60%, -^ 

Mr] 

Ato w 10 ms, A0O ~ 52 rad. 



Axi « 0.4 



(21) 



Evidently, these ranges are extremely large, and we would 
have to incorporate prior restrictions of the parameters to ob- 
tain a slightly more realistic estimate of the parameter uncer- 
tainties [14, 22-24]. However, we merely use it as an illustra- 
tion here. 

Of course, it is well known that particular parameter com- 
binations can potentially be measured much more accurately. 
The canonical example is the chirp mass 



M = M7/3/S 



(22) 



which governs the Newtonian-order PN phase coefficient. In 
the above example. Fisher-matrix calculations in terms of A4 
instead of M (the other parameters remain unchanged) reveal 



AM 



0.32%, 



(23) 



and we shall later formalize the search for such well- 
determined parameters. 
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FIG. 2. The estimated 90% confidence interval around a NS-BH sig- 
nal with component masses 1.35i\f0 and 5Mq, SNR 20 and a BH 
spin of xi = 0.3, as indicated by the (red online) dot. Dashed and 
solid lines are Fisher-matrix estimates in terms of the physical param- 
eters (with parameter differences in terms of A log 77 or linearized as 
Ari/i], respectively). The gray region is obtained by actual wave- 
form overlaps, or equivalently by the PCA introduced in Sec. II. Both 
methods yield indistinguishable results. 



For now, we focus on the other physical parameters and 
show in Fig. 2 the boundaries of the confidence interval, pro- 
jected onto the rj-xi plane according to (17). As we use logi] 
as one parameter, the corresponding variation reads A log 7/ — 
log 77 — log?}, which gives rise to the dashed curve in Fig. 2. 
Since we work only to linear order in the parameters, we may 
as well express A log 77 w A77/77 = (77 — 77)777, which in turn 
leads to the solid ellipse in Fig. 2. Already these two linear 
approximations differ considerably in the range where we ap- 
ply them, indicating that we should not trust the linearization 
for such high parameter variations. 

Instead, the gray shaded region in Fig. 2 can be obtained 
by actual overlap calculations according to (1) and (2), and 
then simply recording the area where the distance between 
target and model waveform satisfies (11). Note that in order 
to compare with the projected Fisher matrix, we optimized the 
inner product over all parameters not shown in Fig. 2 (these 
are to, 4>o and M). This approach, however, entails nontriv- 
ial computational efforts, and a much more efficient method 
is to use the parametrization (15) discussed in Sec. II B that 
embeds the waveform manifold in flat space where the Fisher 
matrix is parameter independent. In that case, calculating the 
waveform difference (4) becomes a simple, yet very accurate 
matrix multiplication. Indeed, when we calculate the data for 
Fig. 2, there is no distinguishable difference between actual 
overlap results and Fisher-matrix estimates in terms of PN co- 
efficients (15), but the computation times differ by several or- 
ders of magnitude, the latter calculation being much faster. 

The details of this fast algorithm are as follows. We calcu- 
late and diagonalize Fjj in terms of the PN coefficients (15). 
Inverting this matrix to estimate the accuracy of these param- 
eters is likely to introduce large numerical errors, because as 
noted in both [17] and [52], large Fisher matrices of this form 
are ill-conditioned and numerical inversion routines cannot be 



Xi Afii (Fisher) A/it (actual) Ai(A/ii)^ (actual) 



1 45300 


0.012 


0.011 


5.56 


2 80 


0.28 


0.27 


5.73 


3 0.84 


2.7 


1.2 


1.17 


4 0.008 


28 


16 


2.08 


5 4 X 10-^ 


390 


1.6 


-10-* 


6 4 X 10-*^ 


10* 


46 


~ 10-* 


7 2 X 10-i« 


2 X 10^ 


4.9 


~ io-« 


8 9 X 10-^* 


8x 10^ 


41 


~ 10-1" 



TABLE I. The principal components of the Fisher matrix that treats 
the PN phase coefficients as free parameters (see Sec. II B). The tar- 
get signal is the same as in Fig. 2. We report the eigenvalues A^ and 
the theoretical spread Afii (Fisher) in the 90% confidence interval, 
assuming all PN coefficients to be free and independent parameters. 
The latter should be contrasted with the actual variation A/ii (actual) 
on the lower dimensional waveform manifold, which in turn affects 
the waveform difference (8) by the amount stated in the last column. 
Note that the numbers shown here depend on the reference frequency 
/o in (13), and we employed /o — 200 Hz. 



trusted. However, diagonalizing the Fisher matrix is numeri- 
cally more stable due to (18) and (19), and we only need to 
accurately calculate the eigenvectors with large eigenvalues. 

Table I reports these eigenvalues A^ that enter the wave- 
form difference through (8). Assuming a maximally allowed 
square difference of 6.25, Eq. (11), we can then simply scale 
the inverse square root of Xi to obtain the theoretical range of 
principal components in the confidence interval [denoted by 
A/ii (Fisher) in Table I]. However, this will be of little value, 
since we cannot easily transform this range back to physical 
parameters, and the actual waveform manifold is only a sub- 
set of the 8-dimensional ellipse we have just calculated, see 
Fig. 1. 

Instead, we densely populate the (physical) parameter space 
around the target signal, transform these coordinates to the 
/ii-space (by a matrix multiplication) and select all points that 
fulfill (11). This is a computationally very cheap procedure 
which allows us to find the actual spread in both physical 
parameters and principal components. Of course, A/i^, re- 
stricted to the physical waveform manifold (labeled by the 
word "actual" in Table I), has to be less or equal than the 
theoretical prediction that assumes all PN coefficients to be 
independent, and indeed, this is what we find in Table I. More- 
over, we conclude from these numbers that only the first four 
principal components contribute significantly to the waveform 
difference, and we can neglect the others for all practical pur- 
poses, as their actual variation is diminished by the small as- 
sociated eigenvalue. 

Let us make a final remark on the power of our approxima- 
tion. Previous uses of PN coefficients as free parameters in the 
Fisher matrix have largely neglected a parameter-dependent 
cutoff frequency in the inner product (2), mainly because the 
considered systems had low enough total mass to place /max 
out of the detector's sensitivity band. For mixed NS-BH bina- 
ries we still may want to neglect the merger and ringdown of 
the signal, but the waveform then has to be terminated in the 
detector band with a consistent cutoff frequency that is at least 



total-mass dependent as in (3). Such additional complications 
do not spoil the accuracy of the estimate, as the following sim- 
ple calculation shows. Assume the waveforms of two systems 
hi, /i2, with total masses Mi < M2. Their distance is based 
on an integral in Fourier space, and due to /max > /max we 
can simply expand 



\hi-h2\\^ = 



/ll - /l2 



\hi\ 



A 



(24) 



where the brackets indicate the integration range, and the sec- 
ond part only contains hi because /12 has been terminated al- 
ready in this frequency range. The first part can accurately 
be estimated with a parameter-independent Fisher matrix that 
incorporates the smaller upper cutoff frequency. The second 
part is proportional to a simple integral over /-^^■^/S'„(/) in 
the respective frequency range. Both contributions have been 
included in the data shown in Fig. 2. 

As evident from Fig. 2, the true confidence interval is con- 
siderably smaller than predicted by "standard" Fisher-matrix 
estimates. From efficiently calculating waveform differences 
for many neighboring points, we can now simply read off the 
range physical parameters that fulfills (11), and we find 



A77 
AM 



< 50%, 



< 40%, 



Axi < 0.25, 
AM 



M ^ ' M 

Expressed in individual masses, we find 



< 0.2%. 



(25) 



2.5M0 < mi < 8.OM0, I.OM0 < TO2 < 2.5M0. (26) 

Hence, at 90% confidence, we would not be able to tell purely 
from a GW observation whether the observed system is com- 
posed of two rather heavy and hardly spinning NSs, or a light 
NS and a significantly spinning BH. The same conclusion was 
recently drawn in a detailed study by Hannam et al. [11]. 



B. Accurately measurable principal components 

Apart from computational benefits, the method of diago- 
nalizing the Fisher matrix, thereby finding uncorrelated pa- 
rameters, enables us to systematically understand which com- 
binations of physical parameters are well measurable and, in 
turn, along which paths GW signals are almost degenerate. 
We consider again our example of a 5M0/I.35M0 BH/NS 
system with a BH spin of 0.3. 

We start with the standard Fisher matrix in terms of the 
physical parameters {log M, log rj, Xi , ioj 4>o}- After project- 
ing out the time and phase shift [see Eq. (17)], the eigenvalues 
we find are separated by 4 orders of magnitude, respectively, 
and the first principal component (with the highest eigenvalue) 
reads 



^1 ex log M + 0.59 log r/ - 0.02x1 . 



(27) 



The spin dependence is rather small, so we neglect it for sim- 
plicity and find that fii ~ e^^ ^ Mrp-^'^ is remarkably sim- 
ilar to the chirp mass (22). It is not surprising, but reassuring 



that the PCA indeed finds the parameter that has already been 
considered as the best-measured quantity as the first principal 
component. Note, however, that the spread of fii in the 90% 
confidence interval is A/ii = Ajji/fli — 0.007%, therefore 
much smaller than the variation in Ai, cf. (23). We can easily 
understand this by imagining the long ellipsoidal shape of the 
estimated interval that extends to very different lengths along 
the principal components. A minimal rotation (such as from 
/xi to M) can dramatically increase the extent of the ellipse 
along the given direction. 

Nevertheless, it is important to keep in mind that the re- 
sults of the above analysis will slightly change with differ- 
ent values of the source parameters, different variants of the 
detector noise curve and other cutoff frequencies. Thus, it 
is likely to be more practical to consider A4 as the approx- 
imately best measured parameter It is still much more ac- 
curately determined than the second principal component /12, 
so we proceed with calculating the Fisher matrix in terms of 
{log A^, 77, XI7 ^0, 0o}- After projection and diagonalization, 
Ai is indeed the dominating contribution to /^i, and /i2 reads 



Ai2 0c 77 + 0.42x1- 



(28) 



We have neglected the small log Ai contribution (that has an 
estimated coefficient of 0.05) in (28) as we regard the chirp 
mass as essentially measured by /ii and we are now inter- 
ested in determining the remaining parameters. As empiri- 
cally observed and discussed in Sec. Ill A, we cannot simply 
determine the individual masses by measuring the two "best" 
parameters very accurately. Even though the chirp mass is 
only mass dependent, the next principal component is a com- 
bination of (symmetric) mass ratio and spin. Thus, as long as 
the spin value is not determined by other means, we cannot 
neglect it. Neglecting it would result in the mass ratio being 
measured incorrectly. 

The second principal component within the 90% confi- 
dence interval is uncertain by A/i2/M2 ~ 1%. which by itself 
is not a large uncertainty. However, to extract the physically 
more interesting parameters r/ and xi and with them the indi- 
vidual masses, we need to consider the third principal compo- 
nent as well, which reads 



^3 ex 77-2.40x1- 



(29) 



We again neglect the small log A4 contribution here (entering 
with a factor —0.02). This third principal component, how- 
ever, is measurable only by An^/pi^ « 200% which severely 
limits our ability to identify the physical parameters. A visu- 
alization of the range of parameters restricted by the spread in 
/i2 and /i3 is shown in Fig. 3. Note that this is entirely consis- 
tent with the standard Fisher-matrix ellipse in Fig. 2 and the 
numbers presented in (21). In particular, we see that due to 
the tilt of /i2 in the 77-xi plane and the fact that 7/ can only 
take values between and 0.25, the allowed spin is actually 
somewhat restricted, whereas we cannot restrict the mass ra- 
tio of the observed system at all, at 90% confidence (and the 
assumptions underlying the detector and waveform model). 

The arguments we have just provided are simple and in- 
structive, but, just as in Sec. Ill A, the results are not very ac- 
curate, and the specific functional forms of /i2 and /ia vary 
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FIG. 3. Illustration of the mass-ratio/spin degeneracy through prin- 
cipal components for the same target signal as shown in Fig. 2. Even 
though fj,2 is accurately measurable (~ 1%), it only restricts the pa- 
rameter space to a line in the r?-xi space. The next principal compo- 
nent, fi3, is less well determined (~ 200%) and restricts the param- 
eter space only weakly. The best measured parameter ^1 determines 
the chirp mass and does not impact this plot. 
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Azi 


Ai2 


A« 


Ai4 


A.5 


Ai6 


A,7 


A,8 




OPN 


IPN 


1.5PN 


2PN 


2.5PN log 


3PN 


3PN log 


3.5PN 


1 


0.98 


0.17 


0.06 


0.02 


-0.03 


0.00 


0.00 


0.00 


2 


-0.18 


0.77 


0.45 


0.17 


-0.36 


-0.07 


-0.10 


-0.07 


3 


0.05 


-0.47 


0.07 


0.17 


-0.65 


-0.20 


-0.46 


-0.25 


4 


0.02 


-0.32 


0.45 


0.25 


-0.27 


0.09 


0.68 


0.30 


5 


0.01 


-0.23 


0.71 


0.07 


0.53 


0.11 


-0.30 


-0.24 


6 


0.00 


-0.05 


0.22 


-0.45 


-0.03 


-0.32 


-0.30 


0.74 


7 


0.00 


0.02 


-0.19 


0.77 


0.15 


0.15 


-0.32 


0.47 


8 


0.00 


0.00 


-0.04 


0.27 


0.25 


-0.90 


0.20 


-0.13 



TABLE II. Contributions of PN phase coefficients to principal com- 
ponents, obtained with a normalization frequency /o — 200 Hz. The 
cutoff frequency employed here is /max ~ 700 Hz (cutoff of the ref- 
erence signal in Fig. 2), although this affects the highest values only 
weakly. 



throughout the parameter space. Again, the better suited 
parametrization in terms of the PN coefficients can poten- 
tially cure both problems to some extent, because we have al- 
ready demonstrated that the Fisher-matrix estimates are much 
more accurate in this higher-dimensional space. In addition 
to that, by using the PN coefficients we automatically impose 
a waveform-adapted functional dependence upon the physical 
parameters that will lead to principal components that are not 
only simple linear combinations of the source parameters. 

Let us stress again that using the PN coefficient (15) as 
free parameters makes the Fisher matrix only weakly varying 
throughout the parameter space, thus we can analyze princi- 
pal components for an entire range of source parameters by 
diagonalizing just one matrix. The transformation matrix Aij 
encodes which PN coefficients contribute most significantly 
to the each principal component, which we illustrate in Ta- 
ble II. We have chosen the cutoff frequency to be that of our 
previous target signal, i.e., (3) with M — (1.35 + 5)Mq, 
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FIG. 4. Contours of the first principal component fii (solid lines), 
displayed in steps of 1000 times the predicted accuracy within a 90% 



confidence interval, A/il 



lOOOA^i (Fisher), cf. Table I. Here 



we assume nonspinning binaries, although the picture is largely in- 
dependent of the spin. Dashed (red online) lines are contours of con- 
stant chirp mass (22) in steps of AM = IMq. 



/max ~ 700 Hz. However, we also tested cutoff frequencies 
between 300 Hz and 2000 Hz (which is even beyond the tidal 
disruption frequency for our example system, but a reasonable 
cutoff for lower mass systems), and the important contribu- 
tions in the first two principal components vary by less than 
10%. 

It is worth pointing out that the numbers in Table II depend 
on the normalization frequency /o chosen in (13), and it is a 
well known ambiguity of any PCA that it changes with scale 
variations of the used variables. We have chosen /o = 200 Hz 
such that Table II gives a good indication of which PN coeffi- 
cients are important in the Advanced LIGO detector band, but 
Aij alone are not invariant quantities. The invariant quantity 
is the waveform difference in the form of Eq. (8), and our aim 
is to illustrate individual contributions to this difference from 
various principal components. Therefore, by visualizing the 
(/o -dependent) values of /ij according to (7), we can immedi- 
ately gauge how strongly the GW signal changes throughout 
the parameter space. 

Fig. 4 shows contours of constant first principal component. 
We plot contours in steps of 1000 times the predicted accuracy 
of /ii in a 90% confidence interval (see Table I), hence we see 
again that /ii is exceptionally well measurable. In addition, 
by simple comparison with constant chirp-mass lines or by 
the fact that the Newtonian phase contribution in ^i is clearly 
dominating, we confirm once more, in a systematic way, that 
the chirp mass is, to a good approximation, the best mea- 
surable parameter in GW inspiral waveforms. However, we 
should keep in mind that the actual best-measurable parame- 
ter is a PN-like expansion series with not only a Al -dependent 
dominant term, but also rj- and xi -dependent higher-order 
corrections. Indeed, Fig. 4 does not change very sensitively 
with varying spin, but we find noticeable deviations of /ii con- 
tours from constant- A^ lines when the spin parameter is close 
to-1. 



In any case, /ii can be very well constrained by GW mea- 
surements, and we use this fact to analyze the other principal 
components in the rj-xi plane only. The other physical param- 
eter, the total mass M, is then determined for each point by in- 
verting ^i(A/, 77, Xi) = const, where we use the value of /ii 
that corresponds to our target signal (see Fig. 2) as the constant 
on the right-hand side. Fig. 5 illustrate the resulting contours 
of /ii, i = 2,3, 4. We find that both fi2 and /.13 are reasonably 
well measurable, i.e., after detecting a signal, we cannot only 
be confident about the value of the chirp mass (under the sim- 
plifying assumptions made here), the associated values of /i2 
and /i3 also restrict the range of plausible source parameters 
to rather narrow bands in the mass-ratio/spin space. However, 
these two bands have very similar structure, and accurately 
identifying the values of 77 and xi individually remains hard. 
This issue is illustrated in Fig. 6, where we overlay the pre- 
dicted confidence intervals of ^2 and 113 around our fiducial 
target signal, and the result we find is entirely consistent with 
the confidence interval depicted in Fig. 2. Note that adding in- 
formation from higher principal components /i,;, i > 4, does 
not add any more constraints as their uncertainty is too large, 
which can be seen for ^4 in the right panel of Fig. 5. 

It is important to realize that the conclusion that three pa- 
rameters are measurable is largely independent of the func- 
tional form of the PN phase coefficients. One might be 
tempted to relate this fact to the three physical parameters we 
started with, but we just showed that three well-constrained 
principal components do not automatically imply that the 
physical parameters can be determined to the same accuracy. 
Even more important is the inverse conclusion. If we had a 
waveform model in the form (13), but with phase coefficients 
that are functions of more parameters (e.g., a second spin or 
tidal parameters of the NS), we would only be able to con- 
strain three parameter combinations (/ii, fi2, Ms)' except when 
the variation of the PN coefficients ^Afc^ through the param- 
eter space of interest is dramatically increased. Of course, 
the functional form of the principal components might be dif- 
ferent, thus Figs. 4-6 may change, but the PCA is performed 
before the phase coefficients have to be specified, so Table II 
and the eigenvalues in Table I are generally valid. 

For convenience and future reference, it is useful to explic- 
itly write down a simplified version of ^2 that includes all spin 
contributions that are known for the relevant terms. Accord- 
ing to Table II, all 8 PN phase coefficients enter /i2 with non- 
vanishing contributions, some of which, however, are much 
smaller than others. We compared the full ^2 expansion with 
simplified versions that included only the one (two, three) 
most dominant A2J contribution(s). The values of such sim- 
plified expressions are different to the full 112 result, but since 
we are interested in variations of ^2 rather than the values 
itself, we only have to make sure that the structure of a sim- 
plified /i2 preserves the one shown in the left panel of Fig. 5. 
(The same argument is underlying our identification of M as 
the approximately best-measurable parameter.) We find in- 
cluding only the IPN and 1.5PN phase terms is not sufficient, 
but by also adding the 2.5PN log term, we find reasonable 
agreement between full and simplified /i2 . 



10 




0.05 0.10 0.15 0.20 0.25 




0.25 




0.25 



FIG. 5. Contours of principal components fii (i — 2, 3, 4), obtained from a PCA in terms of PN ptiase coefficients, see Table II. We show 
slices of constant fii (essentially constant chirp mass). Contours of fi^ and fj,4 are drawn in steps of 2A/ii (Fisher) as given in Table I, i.e., two 
neighboring lines visualize the 90% confidence interval around any point located on the imaginary line centered within this interval. For better 
readability, contours of /i2 are drawn in 5 times bigger steps, i.e., 10A/i2 (Fisher). 
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FIG. 6. Combination of two principal component ranges from Fig. 5 
around the target signal indicated by a (red online) circle. The inter- 
section of both ranges is a good approximation of the 90% confidence 
interval around the target parameters (cf. the actual interval in Fig. 2), 
and it is more accurate than the estimate in Fig. 3. It is still only an 
upper bound because the waveform difference is actually governed 
by the square sum of deviations in /ii (8) and we neglect the effect of 
a mass-dependent cutoff frequency here. 



We therefore conclude 



fi2 w 0.77i/;2 + 0.45^^3 - 0.36V'5 



(30) 



where all PN coefficients are understood as the phase contri- 
bution at /o = 200 Hz. For convenience of the reader, we 
explicitly detail these coefficients below for the more general 
case of two spinning bodies, with spins aligned to the angu- 
lar momentum (recall, the PCA remains unaffected if the NS 
would be spinning, too). In the form used in [41, 42], the PN 
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19x. 
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(32) 



log _ 386457r GStt 
^^ ~ 32256^ ~ 384 



/ 735505 12265 85?? 
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SXa 
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i^xl + xl) 

(33) 



where we used S = {rni — 1112) /M and the spin combinations 



Xs 



Xl +X2 



Xa 



Xl -Xl 



(34) 



It is interesting to note that while 7/'2 is spin-independent, 
^3 and V^g"* contain the leading order and next-to-leading- 
order spin-orbit terms, respectively [35]. The terms cubic in 
the spins are due to the energy flow into the BHs [40] . These 
are less important and not valid for NSs. However, quadratic 
self-spin terms [39] and quadrupole-monopole contributions 
[38] that appear at relative 2PN order (i.e., they are part of 
7/14) affect the overall signal less strongly, as they have no sig- 
nificant contribution to the second principal component /^2- 
The same applies to even higher spin-orbit terms at 3PN or- 
der [36]. We shall verify the importance of particular PN spin 
contributions in the next section properly, but the results are 
already indicated by the form of the principal components. 

We refrain from analyzing /i3 in the same detail. It is 
mainly determined by 7/^5°^, 7/^2 and t\}^^ and is thereby sen- 
sitive to even higher spin corrections. Also the highest or- 
der we consider, 3.5PN {"^t), influences /i3 to considerably 
larger extend than /ii or [i.^- We thus expect that, of the first 
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three principal components, /is will be the most sensitive to 
higher, as yet uncalculated, PN coefficients. This may imply 
that the detectors are indeed sensitive to changes in the values 
of higher order corrections to the PN expansion, even if their 
absolute value is small relative to lower order terms. However, 
once more PN terms have been calculated, they can easily be 
included in our algorithm and the waveform model can be an- 
alyzed for degeneracies with the method presented here. 



IV. MODEL DEPENDENCE AND PARAMETER BL4SES - 
SYSTEMATIC ERROR 
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While the previous section focused on the confusion caused 
by very similar waveforms within one (perfectly known) 
waveform manifold, we shall now turn to systematic errors 
in GW measurements caused by the imperfect knowledge of 
the waveform family itself. Put differently, we shall estimate 
in this section how the recovered source parameters and signal 
strengths are affected when a given signal (the target signal) 
is not necessarily part of the waveform manifold that is em- 
ployed in the search. 

Fortunately, as long as both the target and search wave- 
forms can be expressed in the form of (12) and (13), we can 
still use the linear Fisher-matrix approximation in terms of the 
PN expansion coefficients Qi, Eq. (15), to estimate the effect 
of different waveform families. The only difference to Sec. Ill 
is that now the PN phase coefficients change not only due to 
a change of the physical source parameters, but also due to 

a different functional form -0^°® — i\)\° {Jvl,r\^x\)- This 
means that the following study will be restricted to TaylorF2- 
like waveform representations; however, we are free to modify 
or even drop individual PN contributions to quantify their im- 
portance in a way meaningful for data-analysis applications. 

Our strategy is as follows. Just as outlined in Sec. 11 A, we 
use the Fisher matrix (5) in terms of the PN phase coefficients 
(15), and the transformation to principal components detailed 
in Table 11 is valid independently of the functional form of the 
parameters. Thus, we pick a target signal by fixing the source 
parameters and reference model and transform to the principal 
components as usual. 



Ml 



EA 



n ^3- 



(35) 



We then consider a different search model and transform from 
the associated PN parameters to the same space of principal 
components, such that 



A/Xi = /Zi - /ij 



Ea. 



(36) 



Again, there is no difference to what we did to analyze sta- 
tistical errors, just now dj will be different from Bj for the 
same set of physical parameters. We can, however, vary the 
parameters of the search model to minimize the difference 

min ||A/i|p = mill VAj(A/i,)2. (37) 



FIG. 7. Systematic bias introduced by a nonspinning model search- 
ing for the waveform of a binary with a 1 . 35 Af q nonspinning NS and 
a 5A/q BH with aligned spin as indicated by the horizontal axis. The 
solid (blue online) line shows the bias in the recovered total mass, the 
dashed (red online) line indicates the bias in the symmetric mass ra- 
tio (20). The gray dotted lines shows the bias in the chirp mass (22), 
for which the scale on the right-hand side is valid. 



Note that we effectively also minimize over time and phase 
shifts, but this is implicit in our method through the projec- 
tion of the associated parameters. The difference between tar- 
get parameters and those that minimize (37) are referred to as 
systematic biases, and they indicate to what extent a model- 
based GW search would be confused by the use of an incom- 
plete waveform model. 

As an illustration, let us consider the following scenario. 
The target signal we assume is that for aligned-spin binaries 
including all known spin corrections as detailed in Sec. 11 B. 
This is the waveform model we analyzed in Sec. 111. Fixing 
the masses again to nii — 5Mq, m2 — 1.35Mq, we now ask 
the question how well can the mass parameters be recovered 
if the BH is possibly spinning and the employed search wave- 
form model is that for nonspinning objects ? We easily address 
this question by using standard minimization techniques and 
find the values that minimize (37). Recall, the target param- 
eters define pi and the variable ^i are closed-form functions 
of M and rj (or equivalent parametrizations). As discussed 
above, we do not need to employ sophisticated template bank 
generation algorithms nor calculate direct overlaps between 
any waveforms to answer this question. 

The result of this exercise is shown in Fig. 7. Not surpris- 
ingly, we find that the bias in the chirp mass M is rather small, 
which is due to the fact the the first principal component is 
dominated by this quantity. The second mass parameter, ei- 
ther total mass M or symmetric mass ratio ry, however, can be 
massively biased if the search does not include spin when the 
source is characterized by a considerable spin. We did not re- 
strict the best-matching parameters to physically meaningful 
ranges, and for target spins xi > 0.3 we find that -q exceeds 
its physical range beyond 0.25. We could have anticipated 
this already from the confidence interval shown in Fig. 2 (and 
further interpreted in Figs. 3 and 6), because there the 90% 
confidence interval, extrapolated by eye, would intersect the 
Xi = line at larger, unphysical values of -q. 

Apart from the bias in the recovered parameters, the ac- 
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tual agreement between the best-match waveform and the tar- 
get signal is of interest, as this constitutes an estimate of the 
detection efficiency (i.e., how many signals would be lost in 
the search due to an imperfect agreement between signal and 
search family). As many authors in the GW literature before, 
we shall express the effectualness of the waveform family by 
the fitting factor 



FF — max 



h{e),h{§) 

\hm\\hm\ 



1 — niin 



\h{0)-h{9)\f^ 

2\\h\\^ 



(38) 

where we optimize over the entire waveform family, i.e., over 
all physical parameters 0ph (we added the subscript to distin- 
guish the actual freedom in the waveform manifold from the 
higher dimensional parametrization we are using in this pa- 
per). The right-hand side of (38) can be calculated efficiently 
from the outcome of (37). This equals the fitting factor under 
the assumptions that ||/i(6')|| = \\h{§)\\ = \\h\\ [53-55], which 
is true for our simplified model when we neglect the variable 
cutoff frequency. 

The fitting factor that corresponds to Fig. 7 deviates from 
unity by as much as 15% for highly negative BH spins and less 
than 2% for positive spin values if we allow for unphysical 
symmetric mass ratios. If we restrict 77 < 0.25 then the fitting 
factor drops without bound for increasing spin to an extent 
where we cannot trust our approximation of the inner product 
any more. 

Comparing spinning against nonspinning models was a 
rather extreme case for illustration, and we shall now turn to 
smaller differences between the target model (which we keep 
fixed as the "best" model that includes all known spin terms) 
and the search family. We are particularly interested in the ef- 
fect of various spin contributions to the PN phasing, which we 
will successively drop from the search model to analyze how 
well this "reduced" model can identify the original signal. 

We restrict our study to the case we considered before of 
a I.35M0 NS and a 5Mq BH, but we simulate all BH spins 
—0.95 < Xi < 0.95 in steps of 0.05. For each of these tar- 
get signals we minimize the difference (37) to various search 
models and record the fitting factor as well as the parameter 
biases. Each entry of Table 111 corresponds to the simulated 
spin value that was recovered with the maximal disagreement 
in terms of fitting factor and bias, respectively. 

The search models we consider are as follows: 

no SO tail: Up to 3.5PN nonspinning and spinning contribu- 
tions included (incomplete 3PN and 3.5PN spin correc- 
tions inherited from re-expanded lower-order terms, see 
discussion in [56]), but without the 3PN spin-orbit tail 
contribution derived in [36]; 

3.5/2.5PN: Up to 3.5PN nonspinning and up to 2.5PN spin- 
ning contributions included, i.e., no incomplete spin 
terms considered; 

3.5/2.0PN: Up to 3.5PN nonspinning and up to 2.0PN spin- 
ning contributions included, i.e., next-to-leading-order 
spin-orbit coupling dropped; 

no Xi: Same as 3.5/2.5PN, but without quadratic spin terms 
at 2PN order; 



Model 


1 - FF [%] 


^[%] 


|A,,| 


[%] 


JAxi 


no SO tail 


0.44 


0.20 


46 




0.24 


3.5/2.5 


0.25 


0.06 


28 




0.51 (0.23) 


3.5/2.0 


24 (810) 


0.23 (0.48) 


49 




1.31 (0.63) 


nox? 


0.13 


0.08 


35 




0.32 fO.20) 


3.5/1.5 


13 (450) 


0.18(0.22) 


49 




1.04(0.6.?) 


2.5/2.5 


0.61 (1.02) 


0.50 


74 




1.54(0.60) 



TABLE III. Systematic errors of various models (see text) searching 
for the waveforms that include all known PN terms. The simulated 
target signals are from 5Mq + I.SSMq binaries with the heavier 
object spinning with —0.95 < x^ < 0.95. We report the maximal 
disagreement between the best-fit model waveforms and the target 
signals in terms of the fitting factor (38) and the biases in chirp mass, 
symmetric mass ratio and spin. Italic numbers in parentheses in- 
dicate the values in case we obtain significantly different numbers 
when the search space is restricted to physically meaningful spins 
-l<Xi<l- 



3.5/1.5: Up to 3.5PN nonspinning contributions included plus 
only the leading order spin-orbit coupling at 1 .5PN; 

2.5/2.5: Up to 2.5PN spinning and nonspinning contributions 
included. 

Interestingly, we find from the results in Table III that the 
reduced search models have reasonably high fitting factor 
with the full target waveform if at least the dominant and 
next-to-leading order spin-orbit coupling are included in the 
model. Higher spin-orbit contributions, quadratic self-spin 
terms, quadrupole-monopole interactions and even higher- 
order nonspinning corrections are less important for the de- 
tection of the signal. 

The systematic biases are not as easily interpretable, be- 
cause every search model exhibits almost degenerate regions 
of parameter space with indistinguishable waveforms, similar 
to what we have analyzed in Sec. III. Thus, a template wave- 
form with a much lower or higher parameter bias might agree 
with the target signal almost equally well, but the point we re- 
port as the result of a numerical optimization procedure does 
not include these information. We can, however, compare the 
results in Table III with the statistical uncertainty reported un- 
der the assumption of a perfectly known waveform family, 
Eq. (25), which leads us to the conclusion that the system- 
atic errors reported in Table III are acceptable for the models 
with low fitting factor, except the 2.5/2. 5PN case where non- 
spinning contributions are truncated. 

It is important, however, to point out that neglecting the 
2.5PN spin-orbit coupling leads to a severe loss in FF, and 
searches employing only the leading-order spin corrections 
are prone to miss signals from binaries with considerable 
spin. The numbers in Table III were obtained by allowing 
unbounded values for the spin of the waveform model, and 
indeed, particularly the cases with low fitting factor achieve 
the best agreement with unphysical values of xi < —1- If 
we instead restrict the search parameter space to physically 
meaningful ranges |xi| < 1, we obtain different numbers in 
some cases, given in parentheses in Table III. Note that the al- 
ready badly-performing models then become completely dis- 
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connected to the target waveform space which results in ab- 
surdly high deviation of FF from unity. These numbers are an 
artifact of using the Fisher matrix to estimate large waveform 
mismatches and cannot be trusted. However, the fact that we 
would be unable to detect some spinning systems with such 
models is only emphasized by these results. 

It is interesting to note that we already observed a similar 
effect with nonspinning searches and unphysical values of rj, 
as discussed in connection with Fig. 7. Here we cannot al- 
low for unphysical 77 as some spinning contributions contain 
S = {mi — m2)/M = ^1 — Arj, see Eqs. (31)-(33), which 
has no real solution for ry > 0.25. However, we have just il- 
lustrated that unphysical values of the spin(s) may potentially 
inflate the waveform manifold enough to increase the detec- 
tion efficiency such that signals that are not part of the search 
family have a higher chance of being detected. 



V. CONCLUSIONS 

In this paper, we have considered nonprecessing inspiral 
waveforms of GWs emitted by coalescing NS-BH binaries. 
Such models are essential ingredients for the ongoing efforts 
to directly detect GWs for the first time, and the success and 
astrophysical output of such detections will depend sensitively 
on our understanding of the waveform family employed in the 
search. 

By combining the well-known Fisher-matrix approach with 
a suitable higher-dimensional coordinate choice, we have 
demonstrated that the analysis of degeneracies in the wave- 
form space can be made considerably more accurate than 
previous Fisher-matrix studies of parameter measurabilities, 
while still much faster than full overlap calculations between 
individual waveforms. Even though the high-dimensional 
Fisher matrix may be ill-conditioned, we argued that the 
relevant information about the waveforms can be extracted 
through, instead of inversion, diagonalization of Fisher ma- 
trix. This is because only the eigenvectors with large eigen- 
values affect the waveform considerably. Thus, this procedure 
(which we identify as a PCA) is still well-conditioned, and we 
explicitly presented how we can efficiently find confidence in- 
tervals around a given signal including a parameter-dependent 
cutoff frequency. 

The coordinate choice we employed is based on assuming 
the FN phase expansion coefficients are free parameters [IS- 
IS]. This approach relies on the waveform model being writ- 
ten as a simple amplitude and a complex phase which is a 
sum of purely frequency-dependent functions, each multiplied 
by a single parameter at most. Extending this strategy to ac- 
commodate more complicated functions, such as full inspiral- 
merger-ringdown models or precessing systems, is difficult as 
these models do not obey this simple analytic form. How- 
ever, we restricted ourselves to a regime where the merger and 
ringdown part of the signal do not contribute significantly, and 
recent investigations show that modeling precessing systems 
may be based on a modulation of nonprecessing signals [57- 
59]. In addition, waveform families that are used for detection 
purposes are unlikely to model full precession [46, 60, 61]. 



Thus, even though realistic signals are expected to contain 
some amount of precession, it is worth analyzing nonprecess- 
ing signals first. 

In agreement with previous publications [10, 11, 22, 23], 
we found that the individual masses of the binary's con- 
stituents cannot be well constrained by GW observations 
alone. This is because even though the chirp mass is measur- 
able very accurately, the second mass parameter can be con- 
fused by the presence of spin. Disentangling spin and mass ra- 
tio would require yet another measurement, which is not accu- 
rate enough to place useful bounds on the individual masses. 

With the analysis carried out in this paper, we can now 
rephrase these results in a more formal manner, following the 
results of our PCA. The first, very accurately determined prin- 
cipal component is dominated by the chirp mass (with higher- 
order spin-dependent corrections). The second principal com- 
ponent can be seen as a combination of symmetric mass ratio 
and spin that may somewhat restrict the spin magnitude, but 
does not allow for an unambiguous measurement of either pa- 
rameter. A third principal component is also measurable to 
reasonable accuracy, but it adds little restrictions to the range 
in mass ratio and spin in our case. Higher principal compo- 
nents cannot be well constrained by GW measurements and 
they do not vary enough through the parameter space of inter- 
est to add much information. 

It is important to point out that the explicit form of the prin- 
cipal components is model- and gauge-dependent (in our case, 
the scale freedom is expressed by the normalization frequency 
/o), so the interpretation of the waveform structure in terms 
of principal components reveals no fundamental property of 
the waveform manifold. It is nevertheless a useful concept 
to understand the prospects and limitations of modeled GW 
searches. For instance, we have demonstrated that three pa- 
rameters can be measured accurately, but whether or not these 
lead to astrophysically meaningful statements has to be de- 
termined by the explicit dependence of these well-measured 
parameters on physically interesting quantities. 

This explicit form of principal directions in parameter space 
is, in turn, derived in terms of FN expansion coefficients. 
In our analysis, we have found that higher-order terms also 
have a noticeable influence on the third principal component, 
suggesting that yet undetermined nonspinning and especially 
spinning corrections may influence our ability to measure pa- 
rameters in the future. Among the currently determined FN 
contributions we identified the leading and next-to-leading or- 
der spin-orbit terms as crucial spin corrections that need to 
be included in the waveform model to not change the man- 
ifold drastically. Again, our fitting-factor study of system- 
atic errors was limited to Fourier-domain models of the form 
detailed above. It would be interesting to contrast our re- 
sults with comparisons between various approximants in the 
time and frequency domain. Note however, that time-domain 
models (such as the TaylorTn approximants) can, in princi- 
ple, be transformed to an analytic form in Fourier space as 
well, where the difference between those models and the Tay- 
lorF2 model employed here would lie entirely in undeter- 
mined higher-order phase corrections [62]. Their effect can 
be studied in our framework by allowing the parameter space 
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to be extending beyond 4PN order, which we leave for future 
work. 

As long as such higher-order terms have not been fully cal- 
culated, we need to ensure that the waveform model chosen 
for use in GW searches is capable of detecting signals from 
other, equally plausible models as well. One way of doing 
this is by allowing unphysical source parameters. We have, 
somewhat artificially, compared the full waveform model here 
with reduced search families that were lacking certain con- 
tributions, which we take as a guideline to the situation we 
are actually facing. Namely, that we search for signals in the 
universe (that may or may not be well described by the full 
theory of general relativity) with a restricted, incomplete, PN 
model. We have demonstrated that allowing an unphysical 
spin parameter beyond unit magnitude can indeed reduce the 
systematic difference of waveform families. It remains to be 
tested more thoroughly whether, for the detection problem. 



that reduces ambiguities to a negligible extent. Parameter es- 
timation pipelines, on the other hand, obviously cannot use 
this freedom. However, our algorithm also provides an easy 
way to estimate parameter biases between different waveform 
models with physically meaningful bounds. 
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